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1. INTRODUCTION 

Interest in using solar power plant (SPP) is continuously growing all over the world, despite the 
stochastic nature of solar irradiation, that affects the operation of SPP [1]-[5]. The predominant technology 
all over the world to transform solar energy into electrical energy today is photovoltaics [6]. Solar energy is 
an essential part of modern energetics [7], [8]. The application of solar energy is used in many areas of 
human activity, including stand-alone power systems, to provide electrical energy to small businesses and 
households. A number of scientists’ papers are devoted to the application of SPP in the mineral resources 
sector: i) for stand-alone objects: gas condensate wells in the Arctic [9], stand-alone coal mines [10], [11], 
working sites of oil and gas fields [12]; ii) improving the quality, reliability, and efficiency of the power 
supply [13]-[20]; and iii) principles of creation, management, and economy of energy complexes [21]—[26]. 

The amount of electrical energy that is produced by the photovoltaic panel (PVP) depends on the 
solar irradiation entering its surface, as well as the efficiency of the equipment used to convert and transmit 
energy [27]. Accurate information about the amount of solar irradiation is needed for designing and 
optimizing electrical complexes with PVP. The amount of solar irradiation entering the PVP's surface 
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depends on many factors: geographic location, time of year and day, weather conditions, the quality, and 

cleanness of the PVP, and the orientation and tilt of the PVP. 

High precision and expensive measuring equipment is used to determine the amount of solar 
irradiation in a particular location. For these purposes, a pyranometer, or pyrheliometer, is used. But because 
of financial constraints and/or a lack of necessary equipment, obtaining solar irradiation data might be a 
difficult process, especially in developing states and far-flung regions. Because of this, the role and 
importance of various evaluative and theoretical methods for estimation of the amount of solar irradiation on 
inclined surfaces is growing. 

Hourly [28], [29], daily [30]-[32], monthly [33] models of solar irradiation-depending on the goals 
and tasks of modeling, might be used to design and optimize SPP. For many cases, it is necessary to model 
exactly on hourly time scales, including the time during the full calendar year. This problem is difficult 
because of the stochastic nature of solar irradiation [34]. The solar irradiation that passes through the 
atmosphere is prone to change due to two factors: atmospheric diffusion of air/water/dust molecules and 
atmospheric absorption by ozone/water and carbon dioxide molecules. 

A part of the solar irradiation eventually reaches the earth's surface in the form of diffuse solar 
irradiation [35]. For any specific geographic location for designing and building SSP, knowledge of the 
global horizontal irradiance (GHI) is important, as is the solar irradiance diffused component [36]-[40]. Most 
SPPs have tilted PVP. However, most meteorological stations only have records of GHI, which is why it is 
necessary to convert horizontal solar irradiation into solar irradiation for a tilted surface, taking into account 
its scattered and reflected components [41]. 

The amount of solar irradiation for a tilted surface depends on the clearness index and diffuse 
transmittance index. These coefficients are usually determined empirically for each geographic location. 
There is a high correlation between them, but it 1sn't possible to obtain a unique regression curve for all 
locations or latitudes [42], [43]. The scientific literature about determining the regression dependence on 
hourly time scales has been reviewed by the authors (Table 1 [44]-[55] (see in appendix)) using the 
key words «solar radiation» OR «solar irradiation» AND «diffuse ratio» OR «diffuse transmittance index» 
AND «clearness index» AND «hourly» in the Scopus database. The following conclusions, which are based 
on a critical analysis of the literature review, might be drawn: 

a. The number of works devoted to the determination and refinement of regression dependencies for the 
diffuse transmittance index is large. 

b. Most of the regression dependencies are established based on empirical studies. But there are also 
dependencies that have been established through analytical studies, including dependencies that were 
obtained using data from the NASA POWER database [56], which is freely available. 

c. For Russia and Russian latitudes (4—80?N), the amount of work devoted to the determination of 
regression dependencies for determining diffuse transmittance index, as well as work about the 
determination of the total solar irradiation coming to a tilted PVP, is too small. At the same time, in 
Russia, there are only about two dozen meteorological stations capable of making heliometrical 
observations. 


2. RESEARCH METHOD 

The NASA POWER project application has become widespread among scientists who are involved 
in research on renewable energy sources [57]. NASA POWER (prediction of worldwide energy resources) is 
a tool that provides global information on solar irradiation, wind, and geothermal resources obtained from 
satellites and available for downloading and analyzing. The tool is used for researching, planning, and 
developing energy projects and for evaluating the possibility of using renewable energy sources. The NASA 
POWER database includes long-term estimates of meteorological and solar energy flows that are provided 
by satellite systems. These data meet the needs of scientists and have been shown to be accurate enough 
to provide consistent solar and meteorological data for locations where there is little or no meteorological 
[58]-[60]. 

NASA POWER provides hourly scale information on GHI, air temperature at 2 meters, and wind 
speed at 10 meters. Also, information is provided on the monthly average global irradiance on a tilted surface 
(surface tilt angles: 0°, latitude—15°, latitude, latitude--15?, 90°). Since no information about solar irradiation 
on a tilted surface on an hourly scale is available, a method is proposed for converting hourly GHI via the 
NASA POWER database into hourly global irradiance on a tilted surface. The following parameters 
(Table 2), extracted from NASA POWER, are used as initial data for a model for estimating global irradiance 
on a tilted surface. The proposed method consists of three steps. 
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Table 2. Initial data description from the NASA POWER database 


Parameter 


Description 


All sky surface shortwave 
downward irradiance 
Scales: hourly, monthly 


The total solar irradiation (direct plus diffuse) incident in a horizontal plane on the earth’s surface 
under all sky conditions (GHI). 


2 All sky insolation clearness A coefficient that represents the clarity of the atmosphere is the ratio of the global sky insolation 
index (k;) that is transmitted through the atmosphere to reach the earth’s surface to the average value of the 
Scales: hourly global solar irradiation falling on the top of the atmosphere. 

3 All sky surface albedo The reflection coefficient of the earth’s surface over the entire sky is the ratio of solar irradiation 
Scales: hourly reflected by the earth’s surface to the total amount of incident solar irradiation reaching the 

earth’s surface. 

4 Integrated solar zenith angle The angle between sunbeams and vertical direction. 

Scales: hourly 
5 Top-of-atmosphere The total solar irradiation on the horizontal plane at the top of the atmosphere (extraterrestrial 


shortwave downward 


irradiation). 


irradiance (toa) 
Scales: monthly 
6 All sky surface shortwave 
diffuse irradiance 
Scales: monthly 


The diffuse irradiation on the earth’s surface in a horizontal plane under any sky conditions. 


2.1. Determination of the sun’s position 

Foremost, it is necessary to make calculations to determine the sun’s position on an hourly scale. 
Widely known methods are used to perform these calculations [61]. At twelve o’clock local solar time (LST), 
the sun is at its highest position in the sky. This time is usually different from local time (LT), due to the 
earth’s curved orbit, time zones, and daylight savings time. The meridian of local standard time meridian 
(LSTM) is used as a reference point for a certain time zone, which is determined as (1): 


LSTM = 15?- AT yc (1) 


where ATyrc is the difference in hours between LT and universal time (UTC). 
The equation of time (EoT) takes into account the fact that the earth's orbit and its axial tilt are not 
perfectly round or oriented in the same direction, respectively, and is given as (2): 


EoT = 9.87 sin 2 B — 7.53 cos B — 1.5 sin B (2) 


where B — 360/365(d — 81), d is the number of the day in a year. 
The time correction factor (TC) (in minutes) is used to reflect the change in LST that occurs within a 
particular time zone as longitude changes, and also includes the above EoT. TC is determined by as (3): 


TC = 4(Longitude — LSTM) + EoT (3) 


where Longitude is the local longitude in degrees. 
According to the two previous corrections, LT can be corrected to obtain LST by as (4): 


LST = LT +£ (4) 
60 


The hour angle (HRA) is used to convert the passage of time into a measure of the sun’s position in 
the sky in degrees. At noon, when the sun is in its highest position in the sky, it will always be 0°. The sun 
moves across the sky at an angular velocity of 15° per hour. Before noon, this equates to a negative HRA, 
and after noon, the HRA will be positive. HRA is determined as (5): 


HRA = 15°(LST — 12) (5) 


The inclination of the earth relative to its axis results in a change in the angle of declination, denoted 
by ô, over time. Without this slope, the declination angle would remain at 0°. However, since the earth is 
tilted by 23.45°, the declination can vary plus or minus that amount. On the spring and autumn equinox days, 
the declination returns to 0°. To calculate the angle of declination of the sun, the expression (6) should be 
used. 


5 = 23.45 sin (255) - (d + 284) (6) 


365 
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The angle of the sun’s elevation is the angular measurement of the sun’s position in the sky relative 
to the horizon. The zenith angle is the angle of elevation of the sun relative to the vertical. This angle is taken 
from the NASA POWER database for modeling. Thus, the angle of the sun’s elevation is the zenith angle, as 
indicated as (7): 


a=90°-¢ (7) 


where Cis the solar zenith angle. 

The azimuthal angle from which the sun’s rays reach the earth changes during the day and depends 
on the latitude and time of the year. It is determined by (8). At noon in the Northern Hemisphere, the sun is in 
the south, and in the Southern Hemisphere, it is in the north. 


sin 6-cos Latitude—cos -sin Latitude:cos me) 


arccos ( LST < 12, HRA <0 
cosa 


sin ó:cos Latitude—cos ó:sin Latitude-cos ms) (8) 


@= 


360°— arccos ( LST > 12, HRA > 0 


cosa 


The angle of sunlight incidence on the PVP is determined as (9): 
AOI = arccos(cos £:cos f +sin¢-sinB-cos(@— $)) (9) 
where £ is the PVP tilt angle; ¢ is the azimuth angle of the PVP direction. 


2.2. Estimate of diffuse transmittance index 

The diffuse transmittance index (Ka) is defined as a function of the clearness index (K+). Typically, 
the diffuse transmittance index is determined empirically for every geographic location. In this study, it is 
proposed to define the diffuse transmittance index as a function of the clearness index based on NASA 
POWER statistics data for the observation period 2001-2021. 

The following parameters are retrieved from the NASA POWER database on a monthly 
measurement scale: i) GHI; ii) top of atmosphere irradiance (TOA); and iii) diffuse horizontal irradiance 
(DHD. The values of the clearness index and diffuse transmittance index are found for the monthly average 
values of GHI, TOA, and DHI for each year from 2001 to 2021, according to as (10) and (11): 


GHI 

Ky = FOA (10) 
DHI 

Ka = CHE (11) 


According to the expressions, 252 points are obtained that characterize the dependence of the diffuse 
transmittance on the sky purity index. A regression dependence Kg = f (K¢) builds based on the points 
obtained. 


2.3. Global irradiance on a tilted surface 

The applied model for global irradiance on a tilted surface calculation is based on the isotropic Liu 
and Jordan's model [62]. Liu and Jordan's model is one of the first and simplest models of solar irradiation, 
and this model is recognized as one of the most accurate among the isotropic models of solar irradiation on a 
tilted surface [63]. This model assumes that the intensity of diffuse solar irradiation is evenly distributed over 
the entire sky and is calculated by as (12): 


E = E, + Ej + Ea (12) 
where E; is normal solar irradiation, corrected for the angle of incidence, and is given by (13); E; is reflected 


solar irradiation from the earth and is determined by the expression (14); and Ea is diffuse solar irradiation 
and is determined by as (15). 


E, = DNI : cos( AOI) (13) 
where DNI = s is the direct normal irradiance. 
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E; = 0.5: GHI - p: (1 — cos p) (14) 
Where GHI is the global horizontal irradiance; p is the albedo of a terrestrial surface. 
E4 = 0.5DHI - (1+ cos B) (15) 


Where DHI = Kj, - GHI is the DHI and Kz is diffuse transmittance index (Liu and Jordan model [64 ]). 


3. RESULTS AND DISCUSSION 
The modeling results are presented as tables and figures for two geographic locations: 

Novozapolyarny (Russia) (66.74?N, 79.52°E) and Nefteyugansk (Russia) (61.088°N, 72.6133°E). 

a. Example 1: estimation for Novozapolyarny: the regression dependence definition Kg = f(K;) is 
presented in Figure 1. In Tables 3, 4, and Figure 2, the monthly average values of global irradiance on a 
tilted surface for 2001—2021 are presented by the NASA POWER database and estimated values, 
respectively. Modeling results of global irradiance on a latitude-oriented surface during the year are 
presented in Figure 3. 
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Figure 1. Determination of the diffuse transmittance index as a function of the clearness index calculated by 
meteorological and solar monthly and annual climatology data for 2001-2021 


Table 3. Monthly average global irradiance on a tilted surface for 2001-2021 by NASA POWER, 
kWh/m?/day 


Surface tilt angle 


Mon 0° Latitude-15? Latitude Latitude+15° 90° 
January 0.07 0.66 0.75 0.8 0.8 
February 0.69 1.81 1.96 2 1.98 
March 2.34 4.09 4.25 4.21 4.1 
April 4.49 5.93 5.87 5.61 5.36 
May 5.24 5.54 5.21 4.74 4.41 
June 5.41 5.16 4.59 3.87 3.39 
July 5.13 5.04 4.52 3.84 3.38 
August 3.38 3.72 3.44 3.02 2.72 
September 1.92 2.61 2.53 2.33 2.17 
October 0.68 1.24 1.27 1.25 1.2 
November 0.16 0.75 0.83 0.87 0.86 
December 0.01 0.33 0.39 0.41 0.42 
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Table 4. Estimated monthly average global irradiance on a tilted surface for 2001-2021, kWh/m?/day 


Surface tilt angle 


Month 0? Latitude-15° Latitude — Latitude«15* — 90? 
January 0.07 0.56 0.63 0.67 0.67 
February 0.69 2.10 2.29 2.35 2.33 
March 2.36 4.39 4.54 4.51 4.41 
April 4.48 5.90 5.83 5.64 5.47 
May 5.20 5.51 5.25 4.86 4.59 
June 5.52 5.12 4.56 3.84 3.39 
July 5.13 4.90 4.40 3.73 3.31 
August 3.39 3.62 3.34 2.91 2.63 
September 1.94 2.58 2.49 2.28 2.12 
October 0.69 1.29 1.33 1.30 1.26 
November 0.16 0.73 0.81 0.85 0.85 
December 0.01 0.35 0.41 0.43 0.44 


kWh / m? / day 
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Figure 2. Monthly average global irradiance on a tilted surface for Novozapolyarny 
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Figure 3. Hourly global irradiance on a latitude-oriented surface for Novozapolyarny 
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b. Example 2: estimation for Nefteyugansk: the regression dependence definition Ky = f (K+) is presented 
in Figure 4. In Tables 5, 6, and Figure 5, the monthly average values of global irradiance on a tilted 
surface for 2001—2021 are presented by the NASA POWER database and estimated values, respectively. 
Modeling results of global irradiance on a latitude-oriented surface during the year are presented in 
Figure 6. 
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Figure 4. Determination of the diffuse transmittance index as a function of the clearness index calculated by 
meteorological and solar monthly and annual climatology data for 2001-2021 


Table 5. Monthly average global irradiance on a tilted surface for 2001-2021 by NASA POWER, 
kWh/m?/day 


Surface tilt angle 


Mon 0? Latitude-15? Latitude _ Latitude-15? 90° 
January 0.38 1.08 121 1.27 1.27 
February 12 2:55 2.76 2.83 2.77 
March 2.88 4.38 4.53 4.48 4.25 
April 4.44 5.43 5.3 4.98 4.51 
May 5.08 5.15 4.71 4.11 3.4 
June 5.45 5.14 4.56 3.83 3.01 
July 5.23 5.08 4.55 3.85 3.05 
August 3.75 4.03 3.73 3.27 2.71 
September 2.21 2.8 2.71 2.49 2.18 
October 1.01 1.59 1.63 1.58 1.47 
November 0.39 0.98 1.08 1.12 1.10 
December 0.19 0.63 0.71 0.75 0.75 


Table 6. Estimated monthly average global irradiance on a tilted surface for 2001-2021, kWh/m?/day 


Surface tilt angle 


Month q°_Latitude-15° Latitude Latitude+15° — 90* 
January 0.38 1.50 1.72 1.83 1.84 
February 1.25 2.88 3.15 3.25 3.19 
March 2.90 4.52 4.68 4.64 4.42 
April 4.47 5:35 5.23 4.93 4.49 
May 5.01 4.98 4.55 3.97 3.31 
June 5.39 5.00 4.44 3.74 2.97 
July 5.17 4.92 4.40 3.73 2.99 
August 378 3.91 3.61 3.16 2.63 
September 2.26 2.85 2.76 2.54 2:23 
October 1.02 1.70 1.76 1.72 1.61 
November 0.39 1.12 1.25 1.31 1.29 
December 0.19 1.03 1.19 1.29 1.30 
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Figure 5. Monthly average global irradiance on a tilted surface for Nefteyugansk 
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Figure 6. Hourly global irradiance on a latitude-oriented surface for Nefteyugansk 
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To verify the proposed method, the estimated monthly average values of global irradiance on a tilted 
surface are compared with the monthly average values of global irradiance for 2001-2021 according to the 
climatological data of the NASA POWER (surface tilt angles: 0°, Latitude—15°, Latitude, Latitude+15°, 90°). 
To estimate the effectiveness of the proposed method, the results are compared with the methods of other 
authors (from Table 1), which were designed for the latitude of 60°N and north. Comparisons were held on 
mean bias error (MBE), mean absolute bias error (MABE), and root mean squared error (RMSE) criteria. 
These test methods can be expressed through as (16)-(18): 
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N MBE 


MBE = 45}, i — xi), nMBE = “Z. 100% (16) 


MABE 
x 


RMSE = E» — x;)?, nRMSE ==. 100% (18) 


where N is a number of values; x; is the monthly average global irradiance on a tilted surface for 2001-2021 
by NASA POWER database; x is the average of x;; and y; is the estimated monthly average values of global 
irradiance on a tilted surface for 2001-2021. 

The better the model for estimating hourly solar radiation, the lower the values of MBE, MBE, and 
RMSE [65]. For global irradiance, nMBE and nMABE within +10% and nRMSE<20% indicate good model 
results [66]. Comparison results are presented in Tables 7 and 8: i) model 1 is the model of the proposed 
authors’ method; ii) model 2 is the Reindl et al. [46] model; iii) model 3 is the Berrizbeitia et al. [43] model; 
iv) model 4 is the Bortolini et al. [41] model for an annual scenario; v) model 5 is the Bortolini et al. [41] 
model for a seasonal scenario; and vi) model 6 is the Muneer et al. [55] model. 


MABE = © 5}L,|yi — xil, nMABE = - 100% (17) 


Table 7. Comparison of estimated models for Novozapolyarny 


Parameter Model 1 | Model 2  Model3  Model4  Model5  Model6 
MBE, kWh/m?/day 0.03 0.11 -0.06 0.18 0.16 -0.06 
nMBE, % 1.00 3.96 -2.15 6.53 5.80 -2.15 
MABE, kWh/m?/day 0.09 0.27 0.11 0.23 0.22 0.11 
nMABE, % 3.19 9.61 4.13 8.43 8.06 4.13 
RMSE, kWh/m?/day 0.21 0.85 0.46 1.40 1.24 0.46 
nRMSE, % 7.73 30.65 16.67 50.56 44.89 16.63 


Table 8. Comparison of estimated models for Nefteyugansk 


Parameter Model 1  Model2  Model3  Model4  Model5 Model 6 
MBE, kWh/m?/day 0.09 0.14 0.03 0.25 0.23 0.03 
nMBE, % 3.15 4.73 0.92 8.58 8.06 0.92 
MABE, kWh/m?/day 0.17 0.22 0.15 0.26 0.25 0.15 
nMABE, % 5.80 7.61 5.21 8.98 8.62 5.21 
RMSE, kWh/m?/day 0.71 1.06 0.21 1.92 1.81 0.21 
nRMSE, % 24.39 36.61 7.09 66.48 62.46 7.13 


From the comparison of estimated models for Novozapolyarny, one can see that model 1 is better 
than other models. From the comparison of estimated models for Nefteyugansk, one can see that model 1 is 
worse than model 3 and model 6, but better than model 2, model 4, and model 5. Notably, in model 3 and 
model 6, the diffuse transmittance index was obtained as a result of an assessment of empirically obtained 
data for latitudes close to the latitude of Nefteyugansk. However, Novozapolyarny is located north of 
Nefteyugansk, and in this case, model 1 manifests itself better than model 3 and model 6. 


4. CONCLUSION 

A method for converting hourly GHI to global irradiance on a slanted surface has been developed. 
This method makes use of publicly available data from NASA's POWER database, the well-known solar 
position model, and numerical modeling methods based on an isotropic model of solar irradiation. A 
comparison of the anticipated monthly average sun irradiation values with the monthly average solar 
irradiation values according to NASA POWER climatological data demonstrates good modeling results 
according to the MBE, MABE, and RMSE criteria. 

Furthermore, a method for determining the reliance of the diffuse transmittance index on the 
clearness index was proposed for geographic locations where dependence had not previously been 
empirically demonstrated. This strategy was compared to the models of other authors using the MBE, 
MABE, and RMSE criteria. In modeling, the proposed strategy produces good results. The approach works 
best at latitudes with no experimentally confirmed dependence. However, if there is an empirically verified 
one, it should be used. 
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APPENDIX 
Table 1. Hourly diffuse irradiation models 
No Authors Regression dependency Location Latitude 
1 Orgill and Hollands Ky=1-0.249K, for K, « 0.35 Toronto (Canada) 43.66? N 
[44] Ka = 1.157 — 1.84K, for 0.35 <K,<0.75 
Ka = 0.177 for K; > 0.75 
2 Erbs et al. [45] Ka = 1 — 0.099K, for K, < 0.22 US 31°-42° N 
2 
k = 0.9511 — TIR, + en AE ^ for 022€ K, «08 
—16.638K," + 12.336K, 
Ka = 0.165 for K, > 0.8 
3 Reindl et al. [46] Ka = 1.02 — 0.249K, for K, < 0.3 Europe and North 28°-60° N 
Ka = 1.45 — 1.67K, for 0.3 < K,< 0.78 America 
Ka = 0.147K, for K, = 0.75 
4 Hawlader [47] Ka = 0.915K, for K, € 0.225 Singapore 1.20? N 
Ka = 1.135 — 0.9422K, — 0.3878K,? for 0.225 « K, < 0.775 
Ka = 0.215K, for K, = 0.775 
5 Spencer [48] Ka = a; — b3K, — 0.3878K,? for 0.35 « K, « 0.75 Australia 20°-42° S 
6 Chandrasekaran and Ka = 1.0086 — 0.178K, for K, < 0.24 Madras (India) 13°N 
2 
Kumar [49] Ka = 0.9686 — peas = ee + for 034 « K, « 08 
+10.1862K,° + 8.3733K, 
Ka = 0.197 for K, 2 0.8 
7 Boland etal. [50] Ka = 1/[1 + exp(—5.0033 + 8.6025K;)] Australia (Victoria) 38.09° S 
8 Miguel et al. [51] Ka = 0.995 — 0.081K; for K, € 0.21 Greece, Portugal, 379-44? N 
Ka = 0.724 — PIBER = 8.32K,? + for 0.21 < K, < 076 France, Spain 
+4.967K, 
Ka = 0.18 for K, = 0.76 
9 Oliveira et al. [52] Kyg=1 for K, < 0.17 Sao Paulo 23.54? S 
2 
fuer uu T ake t for 017«K,« 075 
T3AK, + 5.2K; 
Ky = 0.17 for K, = 0.75 
10 Karatasou et al. [53] Ky = 0.9995 — 0.05K, Ta + for 0€ K, < 078 Athens, Greece 37.97? N 
+1.4926K, 
Ka = 0.2 for K, > 0.78 
11 Soares et al. [54] Kg=1 for K, € 0.17 Athens, Greece 37,97? N 
2 
Rupe HAE i ToK * for 017«K,« 075 
0.01K,^ + 3.14K, 
Ky = 0.17 for K, 2 0.75 
12 Berrizbeitia et al. [43] K4 = 0.8636 — 0.9291K, — 0.4623K,? for 13...20°N India, Kingdom of 13°-58°N 
Ka = 1.0815 — 1.8386K, — 0.994K,? for 20...42°N Bahrain, 
Ka = 0.9502 — 1.185K, — 0.8896K,? for 50...58°N State of Kuwait, 
Spain, Portugal, 
United Kingdom 
13 Bortolini et al. [41] Ka = 0.9888 + 0.3950K, — . Europe 37°-59° N 
—3.7003K,? + 2.2905K,? for annualscenario 
Ka = 1.0172 + 0.0158K, — . 
—2.7036K,? + 1.5729K,? for seasonal scenario (summer) 
Ka = 0.9403 + 0.9887K, — i : 
—5.2499K,? + 3.4586K,? for seasonal scenario (winter) 
14 Muneer et al. [55] Ka = 0.95 — 1.185K, — 0.89K,? United Kingdom 51.42°N 
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